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Abstract In this paper we consider the equations of the unsteady viscous, in- 
compressible, and heat conducting magnetohydrodynamic flows in a bounded 
three-dimensional domain with Lipschitz boundary. By an approximation scheme 
and a weak convergence method, the existence of a weak solution to the three- 
dimensional density dependent generalized incompressible magnetohydrody- 
namic equations with large data is obtained. 
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1 Introduction 

The study of the dynamics of biological fluid in the presence of magnetic field 
is very useful in understanding the bio engineering and medical technology. 
The development of magnetic devices for cell separation, targeted transport 
of magnetic particles as drug carriers, magnetic wound or cancer tumor treat- 
ment causing magnetic hyperthermia, reduction of bleeding during surgeries 
or provocation of occlusion of the feeding vessels of cancer tumors and the 
development of magnetic tracers, as well are well-known applications in this 
domain of research [TJ[2] . 

Magnetohydrodynamic flow of a non-Newtonian fluid in a channel of slowly 
varying cross section in the presence of a uniform transverse magnetic field was 
studied in [3]. In the recent past, El-Shehawey et al. [T3] studied an unsteady 
flow of blood as an electrically conducting, incompressible, elastico-viscous 
fluid in the presence of magnetic field through a rigid circular pipe by con- 
sidering the streaming blood as a non-Newtonian fluid in the axial direction 
only. 
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In the present paper, we consider the following system of the three-dimensional 
incompressible magnetohydrodynamic flows 2 , 251134] : 

p t + div(pu) = 0, divu = 0, (1) 

(pu)t + div(pu <g) u) + VP = (V x H) x H + divS(p, 9, D(u)), (2) 

<P t + div(u(#' + P)) = div((u x H) x H + z/H x (V x H) 

+uS(MD(u))+q(p,0,V0)), (3) 

H ( -Vx(uxH) = -Vx(i/VxH), divH = 0, (4) 

where p, u € R 3 , H 6 R 3 and 8 t; S denote the density, the velocity, the 
magnetic field and the temperature, respectively; S is the viscous stress tensor 
depending on the density, the temperature and the symmetric part of the 
velocity gradient D(u), and the thermal flux q is a function of the density, the 
temperature and its gradient; <P is the total energy given by 

<£ = P ( e + i|u| 2 ) + i|H| 2 , *' = p(e+i|u| 2 ) 1 

with the internal energy e(p,6), the kinetic energy ^p|u| 2 , and the magnetic 
energy ^ |H| 2 ; D(u) = Vu+Vu T is the symmetric part of the velocity gradient, 
Vu T is the transpose of the matrix Vu, and I is the 3x3 identity matrix; 
v > is the magnetic diffusivity acting as a magnetic diffusion coefficient of 
the magnetic field. 

Equations ([TJ-Q describe the conservation of mass, momentum, and en- 
ergy, respectively. It is well-known that the electromagnetic fields are governed 
by Maxwell's equations. In magnetohydrodynamics, the displacement current 
can be neglected |2EU22]- As a consequence, Eq. (g]) is called the induction 
equation, and the electric field can be written in terms of the magnetic field 
H and the velocity u, 

E = zA7xH-uxH. 

The basic principles of classical thermodynamics imply that the internal energy 
e and pressure P are interrelated through Maxwell's relationship: 

d p e = \(P- 0d e P), d e e = 8 e Q = c v {6), 
p 2 

where c v {&) denotes the specific heat and Q = Q(9) is a function of 9. Thus, 
we assume that the internal energy e can be decomposed as a sum: 



e(p,9)=P e (p) + Q(9). 



(5) 
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In order to have a clear picture about the admissible structure of these func- 
tions, we assume that S and q are of the form (for p > 0, 9 > 0, D € R 3x3 
symmetric) 

S(p,0,D) = 2p (p, 0,|D| 2 )D, q(p,0,V0) = K O (p,9)V0, (6) 
and behave as 



(7) 



J S(p,0,D) ~ p( P ,0)(e + \D\^-D, rG (l,oo), 
\ q(p,6,V9)~K(p)9 a V6 = ^Vfl« +1 , aeR, 

where e G [0, 1], and there exist positive constants /i, pT, k and 7t such that 

< p < p(p, 0) <JI < +oo, < K < n(p) < ~k < +oo. 

Thus, in particular, for all p > 0, 8 > 0, and D, Be R 3x3 symmetric, 

S(p,0,D) D > M (e+ |D| 2 )^|D| 2 > 0, 
|S(p,0,D)| <7l(e+|D| 2 )^|D|, (8) 
(S(p,0,D)-S(p,0,B))-(D-B)>O, 

and for all V0 G R 3 ( j o < p*), 

f q(p,0,V0)-V0>K0«|V0| 2 = ^|V0^| 2 >O, 
\ |q(p,0,V0)| <K0«|V0|. 

For simplicity, we impose the boundary conditions 

u = 0, H = 0, q-n = on[0,T]xdf2. (10) 

The initial density is supposed to be bounded and the initial total energy is 
integrable, i.e., 

P (o,-) = PoeL 00 (r?), 

p(J^ + d)(o,-) = p ( ] -^f + e )e-L 1 (o), (ii) 

H(0,-)-H o eL 2 (l2), 
and 

< p* < po(z) < P* < +°° a -G- x <E f2, (12) 
< 0* < o (a;) M o.o. (13) 

where p*, p* and 0* are constants. 

The study of long time and large data existence theory for inhomogeneous 
incompressible fluids was investigated in several contributions. For the case 
r = 2 and the viscosity does not depend on |D(u)|, using the concept of renor- 
malized solutions, Lions |30j established a new convergence and continuity 
properties of the density that may vanish at some parts of the domain where 
the viscosity depends on the density. Meanwhile, he got rid of the smallness 
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of the data. For the case that the viscosity depends on the shear rate (r ^ 2), 
Fernandez-Cara et al. ^£5] proved the existence of weak solutions for that 
r > Tf ■ Guillen-Gonzalez [19] (also see [30]) considered the spatially periodic 
setting by using higher differentiability method. Recently, Frehse et al.[T4"] 
established the existence result with non-slip boundary conditions and a vis- 
cosity that depends on both the density and the shear rate for r > ^ . In [T7] , 
Frehse et al. showed the existence result on the full thermodynamic model for 
inhomogeneous incompressible fluids for r > which improved the result 
in [14] . For more results about thermal flows of incompressible homogeneous 
fluids, we refer the reader to [5UHI71IHJE2] - 

If the solution is smooth, multiplying the momentum equation @ by u 
and the induction equation by H, respectively, and summing them together, 
we obtain 

|(ip|u| 2 + i|H| 2 ) + div(^|u| 2 u) + VP • u = divS(p, 9, D(u)) • u 

+ (V x H) x H • u + V x (u x H) • H - V x (>V x H) • H. (14) 

Notice that 

div(^H x (V x H)) = v\V x H| 2 - V x (z/V x H) H, 

and 

div((u x H) x H) = (V x H) x H • u + V x (u x H) • H. 

Subtracting (fT4)) from ([3]), and noticing that divu = 0, we get the internal 
energy equation 

d t {pe) + div(pue) = v\V x H| 2 + S(p, 9, D(u)) : Vu + div(q(p, 9, V6>)). (15) 

It follows from multiplying the continuity equation (TIJ by {pP e {p))' that 

d t {pP e { P )) + dW(puP e (p)) = 0. 

Then subtracting above equation from (|15[) . using ©, we obtain the thermal 
energy equation 

d t { P Q{6))-dw{K {p,9)V9) + dW{ P Q{9)u) 

-S(j>, 9, D(u)) : Vu - v\V x H| 2 = 0, (16) 

where S : Vu denotes the scalar product of two matrices. More precisely, 



S(u) : Vu = 




Recently, there have been much work on magnetohydrodynamics because of 
its physical importance, complexity, and widely application (see [81 l2M2{?ll33j ) . 



Title Suppressed Due to Excessive Length 



5 



Magnetohydrodynamics (MHD) is a combination of the compressible Navier- 
Stokes equations of fluid dynamics and Maxwell's equations of electromag- 
netism. Duvaut and Lions [12] , Sermange and Temam [3 5) obtained some exis- 
tence and long time behavior results for incompressible case. For compressible 
magnetohydrodynamic flows of Newtonian fluids, Ducomet and Feireisl [TTj 
proved the existence of global in time weak solutions to a multi-dimensional 
nonisentropic MHD system for gaseous stars coupled with the Poisson equation 
with all the viscosity coefficients and the pressure depending on temperature 
and density asymptotically, respectively. Hu and Wang [2UJ studied the global 
variational weak solution to the three-dimensional full magnetohydrodynamic 
equations with large data by an approximation scheme and a weak convergence 
method. In |22) . by using the Faedo-Galerkin method and the vanishing viscos- 
ity method, they also studied the existence and large-time behavior of global 
weak solutions for the three-dimensional equations of compressible magneto- 
hydrodynamic isentropic flows ([T])-([3]). They [33J showed that the convergence 
of weak solutions of the compressible MHD system to a weak solution of the 
viscous incompressible MHD system. Jiang, et all. [2~4"1[2"5] obtained that the 
convergence towards the strong solution of the ideal incompressible MHD sys- 
tem in the whole space and periodic domain, respectively. For MHD driven 
by the time periodic external forces, Yan and Li [35] showed that such system 
has the time periodic weak solution. 

The main difficulty of the study of MHD is the presence of the magnetic 
field and its interaction with the hydrodynamic motion in the MHD flow of 
large oscillation. This leads to that many fundamental problems for MHD are 
still open. For example, the global existence of classical solution to the full 
perfect MHD equations with large data in one dimensional case is unsolved. 
But corresponding problem about Navier-Stokes equation was solved in |27) a 
long time ago. In the present paper, we study the existence of weak solutions 
for the density-dependent generalized inhomogeneous incompressible Magne- 
tohydrodynamic flows in a bounded three-dimensional domain with Lipschitz 
boundary. Inspired by the work of [10 , 17, 20 ,31 , we will establish the exis- 
tence of weak solutions for the density-dependent generalized inhomogeneous 
incompressible compressible MHD for any r > 2. 

These equations, and all functions involved in their descriptions as well, are 
considered in (0, T) x i?, where fl C R is an open, connected and bounded 
set with Lipschitz boundary d£2, and T€ (0, oo). 

The paper is organized in the following way: in the next section, by in- 
troducing the appropriate function spaces we provide the precise definition of 
the notion of weak solutions to system (H}-©, ((4]) and (fT6|) . The main result 
of this paper is also stated. Then in Section 3, we first introduce the corre- 
sponding approximation system whose solvability is established in Appendix 
(section 5). We also derive some corresponding uniform estimates. We finish 
the proof of Theorem 1 in section 4 by establishing the strongly convergence 
of {p n }, {u"}, {9 n }, {H™} and {D(u™)}. 
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2 Some notations and main result 

Before giving the definition of the weak solution to the problem ((TJ-©, (0]) and 
(fT6|) with the boundary condition (fT0|) . we first state the following notation of 
relevant Banach spaces of functions defined on a bounded domain Q C R 3 . For 
any p e [1, oo], L P (J?) denotes the Lebesgue spaces with the norm || • Hi/v^n, 

W hp (f2) denotes the Sobolev spaces with the norm || • ||\y 1 '*7f2)' W ' p (i?) 

denotes the closure of C^°(i7) functions in the norm of W 1 ' P (J7). If X is a 
Banach space of scalar functions, then X 3 , X A or X 3x3 , X 4x4 denote the space 
of vector or tensor-valued functions so that each their component belongs to 
X . Further, we use the following notation for the spaces of function with zero 
divergence and their dual (r' = — zj) 

Kdiv^) : = < U 6 Wj^) 3 ;divu = 0}, 

L d . v (L?) denotes the closure of W*'^ (fl) in L q (tt) 3 . The symbols L 9 (0, T; X) 

and C(0, T; X) denote the standard Bochner spaces. We write (a, b) instead of 
f n a(x)b(x)dx whenever ab € L 1 (J7) and use the brackets (a,b) to denote the 
duality pairing for a e X* and b e X. We use C([0, T]; L^, eafe (J?)) to denote 
the space of functions p £ L°°(0, T; L 9 (J?)) satisfying € C([0,T]) for 

all z € L ? . We do not distinguish between function spaces for scalar and vec- 
tor valued functions. Generic constants are denoted by M, their values may 
vary in the same formula or in the same line. 

Definition 1 Assume that S and q are continuous functions of the form © 
satisfying JT])-© with r > 2 and a > — |. The initial data po, u , #o an d 
H satisfy (fTTl) - (fTI?|) . A vector (p, u, 9, H) is said to be a weak solution to the 
problem (HJ-©, (HJ) and ([TB^ of the generalized incompressible MHD equations 
if the following conditions hold: 

• The density 

peC{[0,T\;V(f2)), p t eP 1 (0,T;(W 1,r (i2)r), VgG[l,oo), 

and 

f {p t ,z)-{pu,Vz)dt = Q, (17) 
Jo 

for any z € 1/(0, T; W 1,r (J?)). 

• The temperature non-negative 9 function, the velocity function u and 
the magnetic field H satisfy 

ueL''(0,T;W^. v (r2)), 

5 t (pu)eL'"'(0 ! T;W-^'^)), 
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( P u»eC([0,T];r2) v^4(fl), 
6 e L°°(0,T; L 1 (,!?)), 

e L 2 (0,T;W 1,2 (f2)), 
9 t (p0) e L^O.TjCW 1 ' 8 ^))*) g large enough, 
H £ L 2 (0, T; W 1,a (/2)), H t G L 2 (0, T; W" 1 ' 2 ^)), 

and the following weak formulations hold: 

/ (((pu) t , <p) - (pu ® u, Vp) + (S(p, 0, D(u)),D(p)) 
Jo 

+ [ (H T V(/jH + iv(|H| 2 )-^x))^ = 0,(18) 
Jo 2 

for all^eL f -(0,T;W^. v (^)), 

[ T (((pQ(6)) t ,h) - (pQ(6)u,Vh)-( q (p 7 0,V6),Vh))dt 
Jo 

= f (K|VxH| 2 ,/i) + (S( P ,#,D(u)), D(u)ft)) dt,(19) 
Jo 

for all h G L°°(0,T; W 1,9 (/2)) with g sufficiently large, 

/ ((H i ,i))+y(VxH,Vxi))-(uxH 1 Vxi)p = 0, (20) 
Jo 

for all 6 G L 2 (0,T;W l!2 (!?)). 

• The initial conditions are attained in the following sense 

t Um + ||p(t) - po|lL«(i3) + IK*) ~ u °Hiv 2) + l|HW ~ Ho|l L 2 (r2) = °> V « € t 1 . °°)> 

lim fc) = (p 0o, h), Vh G L°°(r2). 

t — >o+ 

The aim of this paper is to establish the following result. 

Theorem 1 Assume that S and q are continuous junctions of the form (|6]) 
satisfying f^-Q) with r > 2 and a > — g, and t/iere are two positive constants 
c„ and sucft £/ia£ 

< Cu < c„(0) < cr; < +oo. (21) 

Then there exists a weak solution to the problem f2P _ |^P i n the sense of Defi- 
nition 1 with initial data satisfying hll)) - 1iTty) . 
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Note that we consider in (fTS|) only divergenceless test function, so the pressure 
does not appear in the definition of weak solutions. The pressure P (zero 
mean value) can be obtained by comparing two auxiliary Stokes problems 
(homogeneous Dirichlet boundary conditions) 

-Av 1 + VP 1 = div(/m ® u - S(p, 0,D(u)) - (V x H) x H, 
divu 1 = 0, 

and 

- Av 2 + VP 2 = (pu){t) - p u G L 2 ((0, T) x Si) ^ W _1 ' 2 (/2) s , 
divu 2 = 0, 

with taking the test function (p — X(o,t)4 > (X(o,t) denotes the characteristic 
function of (0, t) and G W 'j. (/?)) in (fTS)) . Furthermore, the pressure P 
has the form 

P = P 1 +d t P 2 with P 1 G L r '((0,T) x SI) and P 2 G L°°(0, T; L 2 (i?)). 

The solvability of above Stokes problems can be obtained by a similar proof 
in |^l [T51l37) . Since the presence of dtP 2 , we can not know if P is an integrable 
function on (0, T) x Si. 

3 The uniform estimates 

We take {ipj}JLi as a basis of W '^j v (^) such that (tpi,ipj) = Sij for each 
i,j = 1, ...,oo and ipj G W 'Jj f° r au ii and ^ F n denote the projec- 
tion of L 2 div (S2) onto the linear hull of Let u" G C(0, T; Wj'j. (/?)) 
and H" G C(0,T;Wj'j. (ft)) of the form u n (t,x) = Y^ =1 a^(t)ipj(x) and 
H n (t,x) = Y^j=i c j (t)' t l , j( x ) (note that for r > 2 it is always true that 
W 1 '^ (ft)) ^ W odiv^)) so that the condition divu " = and divH" = 
are automatically fulfilled and (p n , u™, 8 n , H") satisfy 

/ (p?,2>- (p n u n ,Vz)cZt = 0, and p* < < p* in (0,T) X Si, (22) 
Jo 

for any z G L 9 (0, T; W 1,,3 (ft)) with arbitrary o G [l,oo); 

((p"u") t ,^) - (p"u"®u™,V^) + (S(p I1 ,0 n ,D(u")),D(V^)) 

+ / ((H") T V^H" + ^V(|H"| 2 ) • ^)d.x = 0, (23) 
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for alii = 1, . . . , n and a. a. t € [0, T] ; 

7 \{{p n Q(8 n )) u h) + ( P »Q(()>»,V/ 1 )-(q( P " ) f ) Vf),Vf 1 ))<K 



/ (i/(|V x H"| 2 ,/i) + (S(p",r,D(u n )),D(u n )/i))dt,(24) 
Jo 



for all ft G L°°(0,T;W 1,? ) with q sufficiently large; 
The magnetic held function H™ satishes 

(H?, 6) + (V x i/(V x H"), 6) - (V x (u x H"), b) = 0, (25) 

for all b G L 2 (0,T;W 1,2 (J?)); 

0" G L°°(0, T; L 2 (f2)) n 1/(0, T; W 1,p ) with p = min{2, ^^}, 0" > 9, 
in(0,T)xf2; 

The initial data 

p n (0,-) = p , u"(0,-) = r n u O) r(o,-) = o H"(0,-) = H , 

where -T"uo and #q a standard regularization of 9q, fulfill 

r n u — > u Q strongly inL 2 (f2), (26) 
9% — > strongly in L x (/2). (27) 

In fact, using Lemma 3 in section 4, (|22l) is equivalent to 

r2 (p n ,z t ) + (p n u n ,Vz)ds = (p",z)(i 2 ) - (p n ,*)(*i)> (28) 
ti 

for z a smooth function and a. a. t\, ti : < t\ < i 2 < T. 

In the following, we give some uniform estimates with respect to n G N. 
Taking 2 = |u ,l | 2 , t x = and i 2 = < in (gg} we have 

/V, |u"| 2 ) + (p"u", V|u"| 2 )d S = (p™, |u"| 2 )(*) - (p» |r"u | 2 ). (29) 
Jo 

Multiplying the jth equation in (|23[) by Oj, then taking the sum over j = 
1, . . . , n, using (|29l) . and integrating the equality over (0, t), we get 

±{p n , |u"| 2 )(t) + J\s(p n ,9>\ D(u»)),D(u»))da 

+ / / ((H") T Vu"H" + iv(|H n | 2 )-u")^d S = i(p ,|r"u | 2 ).(30) 

«/ vie 



We deal with f|25[) by the same process as in (j30|) . and have 
1, 



2 l|Hn|l L 2 («) 



/ (V x (u n x H"),H n )ds 
Jo 



-jjV* (^VxH"),H")ds + i||r"H ||^ (r2) . (31) 



10 



Weiping Yan, Yong Li 



Direct calculation shows that 

ft ft 



f (V x (u" x H n ),H n )ds =11 (H'fVu"H n + iv(|H"| 2 ) • u n dxds, 
Jo Jo Jn 2 

/ (Vx(i/VxH"),H>= f f v\V x H n \ 2 dxds. 
Jo Jo J n 



So by (f5T|) . we have 



/Vl|VxH"||£ 2 - / ((H"fVu»H"-iv(|H"| 2 ).u»Hx)d s 

= 511^11^,. (32) 

Summing up (|30[) and (13211 . we obtain 

i(p", |u'f)(<) + ^IIH-H^^ + J* HIV x H"||^ 2(r2) + (S(p»,r,D(u»)),D(u"))d S 

= i(p ,|r"u | 2 ) + 5r nH ollL> r (33) 
Then by the first assumption in ((5J), we derive 

(p» |u"| 2 )(i) + HHHI^ + AHIV x Hl^ + (S( P ",r,D(u")),D(u"))d s 

« 

+l£ J* \\B(n n )\\^ (n) d S < (Po, \r n u \ 2 ) + ||r"H ||£ 2(fi) 
<C(/||u || L2(r2) ,||r«H || L2(fi) )<M, (34) 

where M denotes a positive constant depending on the data and maximizes 
all the estimates. 
Note that 

< p* < p n (t, x) < p* < +oo for a.a (t, x) € (0, T) x Q. (35) 
Thus it follows from OU that 



sup ||u n (t)||£ a + sup ||p"|u"| 2 (t)|| Ll + sup ||H"(t)|j 2 <M. (36) 
te[o,T] ^ te[o,T] 1 j te[o,T] ^ 

By Korn's inequality, (|SJ) and ([35]), we have 



< II IF 



+ [ M\ V x Hlf^di + J^W", r , D(u")), D(u»))dt 

+ E £ \\Vu n \\ r Lr{n) d s <M, (37) 

| T ||S(p",r,D(u»))||^, (fi) d s <M. (38) 
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The inequality (|37p implies that 

u" eL r (0,T;W l!r (r2)), H" e L 2 (0, T; H 1 (f2)). (39) 

Using (ESI and (J32J, by Holder inequality and the fact W 1 ' P (L?) ^ L q {(2) 
with p < g < co, we derive 

/ |(p n u n ®u n ,Vu n )|dt< / p*||u"|| 2 ^ \\Vu n \\r»i lo .dt < M,(40) 
Jo 7o L«-i(«) 1 J 

where we require that 1 < pi < r and 

V\ v , „ 

> r =4- pi < < r =► r > 2, (41) 



Pi — 1 r — 1 

which gives one of the restriction for r in our main result. 
It follows from Holder inequality that 



(iFf Vu"Hl Ll(fi) < C(\\U% 2{n) + ||ul Hlw ) < M, 

rni2\ „n\\ s~ r~*(\\ XX n II 2 



l|v(|H«n ■ ui L>) < C{\\n^ W{n) + ||u«|| L2(fi) ) < M. 

Combining above estimates, for 1 < pi < r, we deduce from ([2U]) - ([2~4")) to 

« *<M, / ||ft(p n u n )|f w _ lir , f dt<M. (42) 
(W 'pi-M^))* Jo ( div ( )} 

Let ft = 1 in Using (J35J) and J37|, we have 

sup \\p n Q(O n )\\ L i + sup ||Q(r)|| Ll < M. 

te[o,T] v ; te[o,T] v ' 

By (f2"Tj) and above estimate, 

sup ||p"e n || Ll(r2) + sup ||r|| Ll <M. (43) 
te[o,T] k ; *e[o,T] v ; 

Now we turn to estimate the temperature. Note that 

6 n (t,x) > 8* for a.a. (t,x) G (0,T) x 12. 

Take ft = -{6 n )- x with < A < 1 in (EH). Then by and 



n- A |lL~((o,T )xfi ) < M > ( 44 ) 
IKO^V^H^ dt<M. (45) 



By a contradiction argument, we can easily get 

)^Wv m s cfllHIZ^ + II v ( e»)^- || L , (n) ) if > 0, 



'")^ll L . (n) <c v^|±i<a 
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Thus by these estimates and (|4"4"]) . we have 

rT \wr-^\?^ m dt<M. 

Furthermore, by W 1,2 (/2) ^ L 2 (J?), it holds 

By the standard interpolation of with (|4"5)l . for a > —|, we derive 

jf HHl'^)* < M /or a/Z * e [1, | + a). (47) 
Using ©, (|45j) and gZ]), for 1 < m < we have 



a -J--l-r),dt<M. (46) 



|Ko(p n ,0 n )V0 n | m <ia:cft 

(0,T)xl7 

< « m / |6»™r Q |V6»™rda;dt 
J(o,T)xn 



>(0,T)xQ 

to(q-A-I) m(a + A+l) 

i")-^ — -\ve n \ m {e n ) * — dxdt 

IIV ' 1 M ' L ((0,T)xfi)HV 7 llj^-^-^c^ s \ (0:T)xS - 2) 

< M. (48) 

On the other hand, by ([55)1 and (|4U|) , using the interpolation inequality, Holder's 
inequality and the standard Sobolev imbedding, we have that for some (3 > 1 

T j-T 

".."ni'fl'Mii/ 3 ^ < / p*||u' l||/3 \\a n \\& 
Jo 



\\p n u n Q(n\\$dt < j n P*\\^ n \\i^ {n) \\e n \\^ Hn) dt 

^ cf / T N« n llV-(a)» fln »&-(«)* 



where 



T p 4 3 \ P4 / /.T \ P4 



1 _ 1 1 

P P2 P3 ' 

1/3-^ 1 
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Direct computation shows that 

_ rwj3 2 - rfi _ p 2 /3 _ p 2 (^7/3 - 1) 

P2 ~ mp - /3(r + 1) + ^ ' P3 ~p 2 ~(3' P4 ~ (zo-\)p 2 -vop- 

Note that 1 < p 4 < | + a and a > — ~. We get from (gSJ that 



JO 



y ii u "iiw-J) d N (y ii^h^w* 



where 

p 2 >r, l<P<p 3 <w/3, P<7}, l<m<~. (51) 

Hence we require the restriction that r > 2. 

Finally, by (fM|) . (j3"7| and (|48[). for sufficiently large q, we deduce that 

ii w i nii L >,T;(w-(^) = T iift^nwii^dD).*' 

j o 

= / T sup \(( P n e n ) t (s),h)\ds 

Jo \W\ w ^ (n <i 
< ||fl B WI| L i (fl) < M. (52) 



It follows from ([23) that 
which implies that 



|u n x H n |L j*. <M, 



V x H™ G L 2 (0, T; L 2 (i?)), (53) 



u™ x H" G I/(0,T;L^(I2)). (54) 
Note that ^ > 1 for r > 2. By §U), ([Ml and §5]), we drive 

I^ H 1iAo,t ; (W-(^)).) < M - (55) 

Using Alaoglu-Bourbaki theorem and the uniform estimates (j3"5")l . (|3"6"l) - (f3l))) . 
(|4"3")) - (|47p . (f52")) - (l5"5"l) . for n G N, there exist a subsequences (denoted by itself) 
{p n }, {u™}, {9 n } and {H™} , and (p, u, 9, H) such that 

p" — v p weakly in L p ((0,T) x J?) /or any q G [l,oo) 

and * —weakly in L°°((0,T) x £2), (56) 

< p* < p(t,x) < p* <oo for a.a. (t,x) G (0,T) x 12, (57) 
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fi-Lpt weakly in L Pl (0, T; (W 1 ^ (12))*), (58) 

u"^u weakly in L r (0, T; W 1 '^. (i?)) and L Pl ((0,T) x J?) 3 , 
/or l<pi<r, and * -weakly in L°°(0, T; Lj- V (/?)), (59) 

0"^0 weafc/y in L q ((0,T) x Q) for any q g [1, "t 3 " ), 

o 

n > (9* > /or a.a. (t,x) g (0,T) x 12, (60) 

H ra H weafcZy in L 2 (0, T; W^'j. (/?)) 

and * -weakly in L°°(0, T; L^ iv (L?)), (61) 

VxH"-VxH wea/% m L 2 (0, T; L 2 (J?)), (62) 

H™ — H t weakly in L 2 (0, T; (W 1 ' 2 (X2))*). (63) 
Furthermore, there holds 

p n u n ^pu weakly in L Pl ((0,T) x SI) 3 , 

S(p n ,e n ,B(u n )) ->> S weaA% in L r '((0,T) x J?) 3x3 , 

S = S , S being traceless, (64) 

(6»«)7^(^j? wea^y i n L 2 ((0, T); W 1:2 (J?)) /or 7 € (0, ^y^), (65) 
where pu e L Pl ((0, T)x L?) 3 , S g L r '((0,T)xJ?) 3x3 and 07 g L 2 ((0, T); W 1>2 (J?)). 

4 The strong convergence of {p n }, {u 11 }, {H™} and {Dtu 11 )} 

Before we prove the strong convergence of {p n }, {u n }, {6 n }, {H"} and {D(u™)}, 
we recall two main tools: the Aubin-Lions Lemma (see [36] or [21]) and the 
Div-Curl Lemma (see [HE]). 

Lemma 1 (Lions-Aubin lemma) Let T > 0, p g (l,oo) and {u n }™ =1 be a 
bounded sequence of functions in L p (0,T; X) where X is a Banach space. If 
{ u n}^Li is also bounded in L p {0,T;Y), where Y is compactly imbedded in 
X and {<9tUn}^Li is bounded in L p (0,T; Z) uniformly where X C Z. Then 
{ u n}^Li is relatively compact in {un}^^. 
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For a = (ao, a%, 012, 03) we set 

3 

T>iv t , x a= (a Q ) t + '^2(a i ) Xi and Curl^a = V f ^a - (Vf^a) 2 

i=i 

where V M a = (af , a£ , a^ 2 , a£j. 

Lemma 2 (Div-Curl Lemma) Lt 
Assume that {a"} and {ft™} satisfy 



Lemma 2 (Div-Curl Lemma) Let p, q, I, s £ (l,oo) be such that ^ + ^ = j- 



a n -v a weafc/y m 27((0,T) x i?) 4 , 
6" 6 weafcZy m £ 9 ((0, T) x S2) 4 , 

andDiv t , x a n and Curk, x b n are precompact in W~ 1,s ((0,T)xS2) and W~ 1 ' s ((0,T)x 
J?) 4x4 , respectively. 
Then 

a" • 6" a - b weafcZy in L l ((0,T) x S2), 
where ■ represents the scalar product in R 4 . 
By ([56")) - (l5Tl) and Lions- Aubin Lemma 1, we have 

p n — >p strongly in C([0, T}; (W 1 '^)*). 

We take a™ = (p n , p n t%, p n u%, p n uf) and b" = (^,0,0,0), i = {1,2,3}. It 
follows from (1551)- (15711 that 



a™ — 1 (p, pui,pH2,pus) weakly in L 9 ((0,T) x Q) 4 V q £ [l,Pi], 

b" 0,0,0) weafc/y m L Pl ((0, T) x SI) 4 . 

Furthermore, we notice that L r ((0,T) x L?) 3x3 <-><-> W~ 1,r ((0,T) x SI) and 
Div M a" = (p") t + div(p"u") = 0, 

( Vu" \ 

Curl t 2; b™ = „^ T I , (o denotes zero 3x3 matrix) 



-(Vu" 

and Vu™ is bounded in L r ((0,T) x J?) 3x3 . By Div-Curl Lemma 2, we have 

p n vZ^pUi weakly in L q ((0,T) X Si) Voe[l,^-]. (66) 

It follows from (|58p and (1551) that p and u satisfy (1T71) . Take the test function 

1 P1 

of the form X(ti,t 2 )' l i h £ W 'pi- 1 (J7) in (flTl) . Then partial integration with 

1 1 P1 

respect to time and the density of L ((0,T);W implies that p £ 

C([0,T];L~ afc (/?)),i.e., 

lim (p(t),fc) = (p(i ),/i) Vt € [0,T]. (67) 

t — >t 
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Using the concept of renormalized solutions to the equation (|2"2"j) and following 
the method in [31] . we have 

p n — ► p strongly in C([0, T]; L q (S?)) Vq G [1, oo) and a.e. in (0,t) x Q, (68) 
and 

f hm + ||p(t) - po|lL«(n) = G [1, oo). (69) 
In the following, we prove the convergence of {u™}. By (pL2")) and (IfJBl . we have 
9 t (p n u") - d t (pu) weakly in L r '(0, T; W^'(/2)). (70) 
Furthermore, it follows from (|36|) . ((59)) and f|68[) that 

Vp^u™ loeafcZy in L 2 ((0,T) x ttf. (71) 

Introduce the Helmholtz decomposition 

u = H[u}+H ± [u], 
H ± [u}=V(f>, H[u]=curl<p, 

where <fi is given by the solution to the Neumann problem 

Acf> = divu x € SI, 

^ = o x e 9i7, [ 4>dx = {), 

on Jn 

and <p satisfies the following elliptic problem 

curlH^u] = Vu = w, x e Q, 
div% x [u] = 0, xei7, 
"H^u] • n = x e df2. 

Then by the Aubin-Lions Lemma 1, (g2J|, <f7T|) and ([70]). we have 

(p"u") div — > (pu) div siron^ *n C(0,T;W^' (/?)), (72) 
which together with (|59|) implies 

y/p^u n (t) — > y/pu(t) strongly in L 2 ((0,T) x H) 3 for almost all t £ [0,T]. (73) 
This together with (g^J) and ([55]). ([57 )1 - ([59" ]) shows that 

u" — > u strongly in L 9 ((0,T) x J?) 3 for all q € [1 ,Pi) and a.e. in (0,T) x J?, 
and for r > 2 

p n u n <g> u" ->> pu ® u weafc/y in L r '(0, T; L r ' (J?)). (75) 
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Now we give the strong convergence properties of {#"}. We take a™ = (p n Q(8 n ), R" , R% , R§) 
withi?" := p n 6 n u n + n(p n ,9 n )V0 n , and b™ = ((0™)^, 0, 0, 0) with /3 G (0,^) 
rather small. Then by (gHJ), ® , (|S0 ]l and ([68 ]) . we observe that for ± + ± = 1, 

a™ weakly in L a ((0,T) x Q), 

b n ^(¥, 0,0,0) weafcZy m L k ((0,T) x ft). 

For 1 < si < 2, we have 

Div M a" = d t (p n Q(6 n )) + diviT = u\V x H'f 

+S(p n ,6 n ,T>(u n )) -D(u") c L 1 ((0, T) x 12) m-m- W~ Ml ((0, T) x fi), 

and 

Curl t >" = (_ (v( °„ )7) t V( ^" )7 ) C L 2 ((0,T) x X2) 4x4 ^ W^ 2 ((0,T) x ^) 4 * 4 . 

Using the Div-Curl Lemma 2, for some <7i, <T2 > 0, we have 

p n {9 n ) 1+1 ^ pOW weakly in L 1+CTl ((0, T) x Q), 

p(6 n ) 1+ "< p9Wy weakly in L 1+CT2 ((0, T) x Q). (76) 

Then by Minty's method and ([TBI , we get 

F= 7 a.e. m (0,T) x Q. (77) 

It follows from ([76 | -([77 )l that 

p^T6»" — ► p^T6» strongly in L 7+1 ((0,T) x f?), 

which together with (H7|) . ([57| and (|60l) implies 

0" — >9 strongly in L 9 ((0, T) x Q) M q a). (78) 

By dS3), (JS5J, {BS]), CU) and (JED, we have 

^"0™ — ■» p(9 strongly in L 9 ((0,T) x fl) Vq G [1, - + a), (79) 

o 

p"fu" — -> pflu strongly in L 1 ((0,T) x 12), 

(6»«)7^g|7 ujgafc/y in L 2 ((0,T);W 1,2 (fl)) /or 7 G (0, ^y^), 
which implies that 

V(6»") ii ^ ±i ->> Vtf 2 ^ 1 weakly in L 2 ((0, T); L 2 (J?)). (80) 
It follows from (JUJ) and (JTSJ) that 

p n Q(8 n ) — > P Q{8) strongly in L q ((0,T)xQ) Vge[l,-+a). (81) 
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Note that 

q.(p n ,e n , ve n ) = K (p n ,d n we n = — - — (e n ) 1± ^K O (p n ,e n w(0 n ) s ^ ±1 .(82) 

a — A + 1 

Take A > sufficiently small and q > 2 satisfying ( a+x + 1 '>i = | + a — A. Then 
we have 



T 



o L ^ JO JO 



W n ) 1± ^K Q {p n ,e n )\\\ qin ,dt<C / (fl") 121 ^^ < M. (83) 



Using Vitali's theorem, (|68|l. ((78)) and (|83)l. we derive 

(0 n ) i± T :SL K O (p n ,0 n ) — -> (^) il T =2 KoO',^) strongly in L 2 (0, T; L 2 (L?)). 
So by (@B|), (EOD-dMl), we get 

q(p n ,0 n ,V0 n ) — -> q(p,0,V0) weakly in L*((0,T)xf2) Vz G (1, *±*EL\ 

4 + 3a 

It follows from flU]) and (J53J) that 

H" — ► H strongly in L 2 ((0,T) x ft). (84) 

Furthermore, by (|59l) and ((84)) . we obtain 

u"xH"^uxH weakly in L 2 (0, T; L 2 (ft)). (85) 

Finally, we study the convergence property of {D(u")} by applying the Minty 
method (standard monotone operator technique). The following integration 
by parts formula (the case r > 2) can be obtained by a small modification of 
Lemma 4.1 in [T7]. For reader's convenience, we give the proof. 
We define 



1 f h If 

* z)(t,x) := — / z(t + T,x)dr, (u h * z) (t, x) := — / z(t + T,x)dr, 
Jo " J-h 

. _ z(t + h,x) — z(t,x) h _ z(t,x) — z(t — h,x) 
h ' Z '~ h ' 

where z denotes any locally integrable function and h > 0. 
Then we have the following relationships: 

(u+*z) t =D h (z), (uj^*z) t = D- h (z), (86) 
- T f(r)D h g(r) = t D~ h f(r)g(r)dr + \ f 



h 



to+h 

f(r - h)g( T )dr 



ti+h 

f(r - h)g(r)dr.(87) 

H 
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Lemma 3 Assume that r > 2, 



p G L°°((0,T) x[2)n C{[Q,T];L"{Q)), Vg G [l.oo), 
u G L r (0,T; W^ div (f2)) n £°°(0, T; L 2 (Q)% 



d t (pu)G(L r (Q,T;W^ div (n)) 



and the coupled (p, u) is a weak solution to p t + div(pu) — 0, which means that 
for allze £ r (0,T; W^* (!})), z t G L 1+a (0, T; L 1+a (f2)) andV t a ,h G [0,T], 

ti 

(p(s),z t {s)) + (p( a )«(s),V«(«))da = (p(ti),z(ii)) - (p(t ),z(t )). (88) 

to 

TTien the following formula holds 
r*i 

(9 t (pu),u)-(/5u®u,Vi0di = if(*i)-#(io), V to, ti G [0, T], (89) 

to 

where K(t) is defined by 

K{t) := \{p{t), \u{t)\ 2 ) = 1 -J^p{t 1 x)\u{t,x)\ 2 dx = \\\{Jpu){t)\\l 



Proof To prove this Lemma, we need 



59i T 

r^ Pl = -, r > 2. (90) 



pi — 1 r — 1 

Define 

rti 

L h := I ((pu)t, w+ * ui h *u)dt, 

Jt 

where < t < t x < T and Vft G (0, min{T - i 1( t }). 
Integration by parts on (to,fi), we have 

L h = - (pu,D h (oj- * u))dt + {{pu)(h), (cj+ * uj- * u)(ti)) 

- ((pu)(* ), (wji * * u)(t ))dt. 



Thus to prove (|8§)). using (|87p . it is equivalent to show that 

'to 



lim L h = [ 1 (pu, iv|u| 2 )dt + [^(tj) - A-(to)] 
fc->-0+ J t , 2 



r*i 

lim / (pu,D' l (^ *u))dt + 2[Jf(ti) -If (to)] 



fi->0+ 



to 



't 



lim / (D~ h (pu),uZ *u)dt. (91) 

^^0+ Jtn 
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By and D h p = — div(cj h * (pu)), we derive 
f ti 

/ (D- h 0ni) >W fc *u)dt 
= (pD- h u,uj' h *u)dt+ / (L>" /l p)u(--^),w,;*u)dt 



to 



(p. ^ * u l*) rfi+ / H *(p u )> v ( u (--^)-K *u))dt. 



Finally, due to ((90)) . we can take z = \\w h * u| 2 in ((88]) . Then we conclude 
that 



rti 

/ {D- h {pn),w^*n)dt 



= (p(*l), \K * ll| 2 (fl)) - (p(t ), l|u£ * U| 2 (t )) - jT '(PU, l| W fc * u|f )dt 

+ / * (p u ), v(u(- - • ( w h * u ))^- 

J t 

Thus taking limit h — > in the above inequality, and recalling (|91[) . we can 
obtain (l89|) for almost to and t\ in (0, T). 



Using jBT]), d64|), dZQJ), dZSJ, (|62 )l -(|63 l) and (|85J), we take the limit n — > oo in 
(1231) and © and get 



T 

({pu) t , ij>) - (pu <g> u, V^) + (S, D(V0) 



+ ^((H) T V^H+iv(|H| 2 )-^)dx^ = 0, VVeW^. v (/2), (92) 

y (H t , 6) - (V x (u x H), 6) + (V x (i/V x H), 6)«ft = V6 e Wj*J iv (/2). (93) 

Let X(t .ti)( s ) be the characteristic function of the interval (to,t\) and € 
W 'V (J?). We take the test function ip{s,x) = X(t .ti)( s )'/ , ( a -) m (|89p . then 
performing partial integration with respect to time in the first term, for V</> £ 
Wjj. v (/2),weget 

(p(f)u(i),<£) - (/OoUo,<^) 

- jf (pu ® u, V0) - (S, D(0)) - (H) T V0H - iv(|H| 2 ) • <f>dt, 



which implies 



lim >(t)u(t) - p u , 0) = e L^j (ft). (94) 
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By (jT5|) . letting n — > oo in (|3"0")) . we have 

(p,|u| 2 )(t) < (po,|u | 2 ) /or a.a. te [0,T]. (95) 
It follows from (JBTJ), (JM^-dlS]) that 

t lim + || y/rf)(u{t) ~ "o)||^ ( , 2) = t hm + ((p, |u| 2 )(<) - 2((pu)(t),u ) + |u | 2 )) - 0. 
This together with implies that 

lim>,|u| 2 H(p ,|u | 2 ). 

t — >o+ 

Then by Lemma 3, there exists T a G (0, T] such that 

° (d t ( P u),u) - (pu (8 u, u)dt = |u| 2 )(T ff ) - i( Po , |u | 2 ), (96) 

where T a — > T as a — > 0. 

Taking ip = X(t 0) ti) u an d & = X(t ,ti)H in (j9"2l and (1§3"1) . by we derive 

|u| 2 )(T CT ) + £ ((S, D(u)) + ^((H) T VuH + iv(|H| 2 ) • ujds) eft = i(po, |u | 2 ), (97) 

and 

\^{T a )fv {Q) + ^H|V x H||£ S(fl) - jf ((H) r VuH- iv(|H| 2 ) ■ u)cfe)d* 



^llHoll^^. (98) 



Summing up ((97)) and (|98| . we obtain 



2 



'L («) 



rff 



i((p, |u| 2 )(T CT ) + IIH^JH^^) + ^ ((S, D(u)) + HIV x H 

= ^((po,|uo| 2 ) + ||H || 2 - 3(n) ). (99) 
By ®, for B e L r ((0,T) x f2) 3x3 and B = B T , we have 

/ (S(p n , 8" ,D(u ra )) - S(p n , 0" , B), D(u n ) - B)dt > 0. 
Jo 

Using (|33[) . we can rewrite the above inequality as 

o < \{{ PQ , |r-u | 2 ) - or, |u-| 2 )(r ff )) + ^(IIH"^)!!^^ - im 2 . a(J3) ) 

+1/ / " ||V x H n || 2 di- / CT (S(p",r,D(u")),B)d< 
Jo Jo 

-/ (S(p n ,0",B),D(u n ) _B)dt. (100) 
Jo 
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Let n — ► oo in ([TDD]). By ©, (pp])-([27]). (JM1), dSHJ) , CZEJ, (HI and the 

Lebesgue dominated convergence theorem, we derive 



i(( P0 , |uo| 2 ) - (p, |u| 2 ) ( r CT )) + i(l|H(T ff )||^ (r2) - ||H |£ V) ) + „ jf ' ||V x H|||, (fl) 

- / "(gJ.Bjtft- / "(S(p,0,B),D(u)-B)tifc>O, 
Jo Jo 

which together with (IPPj) and letting cr — s- + implies 

/ (S- S(/9,6»,B),D(u) - B)di > VB e L r ((0,T) x r?) 3x3 , B = B T . 
Jo 

Taking B = D(u) ± tC with any symmetric matrix C E L r ((0,T) x fi) 3x3 
and t > 0, then using Minty's method, we have 

/ T (S-S( P ,tf,D(u)),C)dt = 0, 
Jo 

which means that 

S(<, x) = SO, 61, D(u))(t, x) /or o.o. (t, x) € (0, T) x Q. (101) 

This together with ([53|1. (|D2j) shows that {15]) is satisfied for p, u and H. 

Next we show that jTSJ holds. It follows from (gD]), d73J), O and 

(fTUTj) that 

lim / T (SO", r,D(u")), D(u")) - SO, 9, D(u)),D(u))) df 

= 1 lim lim (0, |u| 2 )(T ff ) - 0", |u"| 2 )(T CT ) - Oo, |u | 2 ) + (p, |^"u„| 2 )) = 0. 

By ©, ([Ml), dSHl), (EH, C3) and §H]), using the Lebesgue dominated conver- 
gence theorem, we have 

lim / (SO",6»",D(u))-SO,6 , ,D(u)),D(u n -u))rfi = 0, 
lim / (SO,0,D(u)),D(u"-u))di = O, 

n — loo],, 

lim / (S(p n ,0 n , D(u")) - SO, 0, D(u)), D(u)) dt = 0. 

n >oo _y 

Above four convergence results show that 

lim / (SO n ,6»",D(u n ))-SO",6i",D(u)),D(u"-u))rfi = 0, 

n >co J Q 
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which combing with shows that 

< (S(//\(9 n ,D(u")) - S(p",<9'\D(u))) ■ D(u" - u) — >0 strongly in L 1 ((0, T) x Q). 
Consequently, for any h € L°°((0,T) x fl), we get 

0= lim / S(p n ,6 n ,r>(u n ))-S{p n ,e r \T>{u)))-T>(u n -u)hdxdt 

>00 J(o,T)xa 

= lim / S(p n ,e n ,B(u n )) ■ B(u n )hdxdt 

n ^°° J(Q,T)y.n 

+ lim / S(/9™,6» n ,D(u)) -T>(u)hdxdt 

n ^°° J(0,T)X(2 

- lim / S(p n ,6 n ,T>(u))-T>{u n )hdxdt 

- lim / S(p n ,6 n ,D(u n ))-T>(u)hdxdt, 

n— *-o°y(o,T)xfi 

which together with d64}, (EH), (EH, dZHJ) and ([84]), the Lebesgue dominated 
convergence theorem implies 

lim / S(p n ,6 n ,T>(u n )) -T>(u n )hdxdt 

n ^°° J(o,T)xn 

= lim / S(p,0,D(u)) ■ ~D(u)hdxdt. 

Thus it follows from the above convergence properties, (|8ip and letting n — > 
oo in d2U) that ([11]) holds. 

5 Appendix-Existence of solutions for the approximation system 

In this section, we establish the existence of solutions to the approximation 
problem ((T])-©, (IT51) and Q for any fixed n £ N. We introduce a two-level 
approximation depending on the parameters e > and k € N, then study the 
behavior of the system for e — > and fc — > oo. 

Let {wjjj^i and be a smooth basis of W ' (J?) orthonormal in 

the space L 2 (i7) such that (w^ojj) = 5^- and (ra7i,n7j) = S^j. For e > 
and fixed k <E N we seek approximation solutions (/O fc ' e , u fc ' e , 9 k,t , H fc,e ) where 
(u fc ' e , 6» fe - e , H fc ' e ) has the form 

j=i i =1 i =1 

and solve the following approximation system (for all 1 < j < k and a.e. 

ie(0,T)) 



p t fc,e + div(/' e u fe ' e ) - eA/- 6 = 0, on [0, T] x [2, 



(102) 
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(p^lu^,^) + (^VV'^i) + (S fc>e ,D(^)) - e(V/' e , [Vu^]^) 

HfVxtf-'jxH 1 ^,), (103) 



(p*. e -Q(0 fc . e ), Wi ) + ( J fi*Q{6 k <*)xL k >', Uj ) + (q( Pk , e ,6 k >*,ve k >% v Wi )) 



= (|V x H fc > e | 2 , W j) + (S M ,D(u fc '>,), (104) 



(4 Hfe,e ) + i/(curlH fc ' e , curlnjj) - ((u fc < £ x H*' e , curler) = 0, (105) 

with 

p fc - £ -n = on [0,T], p* < p fc - £ < p* , p fc ' e (0,-)=p in (0,T) x Q, 
u fe < £ (0, •) = u k ' £ = r fe u , fe < £ (O, •) = 9^ = r k 9 , H fc ' £ (0, •) = H fe ' e = r fc H , 

where 

S fc , £ :=8{ffi e ,6!& x) I> k ' e ) with ^:=max{^,0}. 

The existence of solutions for the approximation system (|102j) - (ll05p can be 
proved by modified Faedo-Galerkin method. We first give the solvable of in- 
duction equation (|105|) . 

Lemma 4 Assume that the initial data TS^'^O) G Y* and given a veloc- 
ity field u k ' t G C([0 1 T},X k ). The system has a solution B^' e (x,t) € 
C^QO, w]; W 1,2 ). Moreover, the operator u k > e — > H k ' e (u k > e ) maps bounded sets 
in C([0, T], X k ) into bounded subsets of z, and the solution operator is con- 
tinuous operator. 

Proof Define 

k k 
X k := {u*'V e 0M) = J24 € Y fe := {H fc » e |H fe ' e (z, t) = £ cf (t)^(x)}. 

3=1 3=1 

Note that H fc ' c (i,x) e Y fc . We can write 

k 
3=1 

where the coefficients c k ' e (t) are required to solve the system of ordinary dif- 
ferential equations 



dC 3 

dt 

\3\<k 



£^(t)cf = 0, b1<fc, (106) 
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where 

Aj. k (t) = v(Vvji, Vruj) - ((u fc - e • - {{wi ■ V)u k ' e ,TUj) 

+ ((divu fc ' e )tn i , vjj ) . 

For given initial data H fe,e (0) <E Y fc , the system (| 106[) has a unique solution 
<y e G C 1 ((0,T); Y fc ) for some T < T. Multiplying both sides (fT06j) by c k '\ 
summing over j, integrating by parts, we have 

J t \\K k ' e Wb im + MVH M IIL- (1J) = -2(H fe - • Vu*-,H*-). 
Note that by Young inequality, 

- (H„ • Vu fe,e , H fc,e ) < ^\\n k %* w + ^llVu^-Hji.^, (107) 
Using the Poincare inequality, 

IIvh*!^ > c\\n k ^ {oy (108) 

By ([TOT jl -lfTng ]) . we derive 

|l|H^||i 2(r2) + cH|H-||^ (f2) <A||vu-|| L2(r?) . 

Thus, for t G [0,T], we get 

m k ' € \\h w < m k ' € (0)\\h {a f- Cvt + ^ [ e- c ^||Vu^|| L2(f2) d s .(109) 
This implies that for t G [0,T], 

HH^IIL^) * HH^(0)||i V2) e-^ + ^ jT e-^-)|| Vu^|| LV2)ds . 

Due to {zuj(x)}j^ =l be an orthonormal basis of W 1,2 (J?), so we have |c^' e (t)| = 

||H fc ' e (t)B 2 , from which we conclude that T = T. 

Define the ball Br of radius R and the map 77 : Br — > Br such that 
77(H fe ' e (0)) = H fc ' e (T), where the radius 7? such that 



7? > 



1 - e 



-CVT 



Follows |31) . we can prove the map 77 is continuous. Hence, it has a fixed 
point. Moreover, from (I109p . we know that the solution operator u k ' e — > 
H fe ' e (u fc ' e ) maps bounded sets in C([0,T],X fe ) into bounded subsets of the 
set Y fe . Then, as done in [32], the solution operator u fc,c — ► H fc ' e (u fc ' e ) is a 
continuity operator. This completes the proof. 
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The rest process of proof is similar to Proposition 7.2. in fTJ] or Lemma 3.2. 
in |22j . Thus, combining with Lemma 4, we obtain the following result: 

Lemma 5 Under the assumption in Theorem 1, for fixed k 6 N and e > 0, 
the approximation problem ilQ2\) - !!TUS\) has a solution (p fc ' e , u k,e ,0 k,e , i2*" ,c ) on 
(0, T)xfi for any fixed T > 0. 



5.1 Limit e — > 

First, we summarize the estimates available for (1102|) - (|105[) for e > and k £ N 
fixed. Then the behavior of relevant solutions will be studied as e — > 0. 
Multiplying (fTU2l by p k > e leads to 

j n^iiLv) + 2e f "^"IIl'w* * Mi,-™- ( n °) 

It follows from (|12j) and a weak maximum (minimum) principle that 

P* < < P*. (HI) 
Taking the L 2 scalar product of (|102[) with a smooth z leads to the equation 

(p k >% z) - (p k <*u k <*,Vz) + e(Vp fc < e , Vz) = 0. (112) 

Multiplying the jth equation in (|103l) by a^' £ , then taking the sum over j — 

1, . . .,n, using (jl 12[) with z = J — and integrating the equality over (0,t), 
we have 

W" k ' e \\h im + ll^» M llL»(fl) + 4 jf (S^D(u^))*. 

+ / / (H fc < e ) T Vu fc < e H fc < £ + iv(|H fe < £ | 2 )u fe ' £ dxds 

< 2||V^r"u ||£ 2(fi) . (113) 

Multiplying the jth equation in (|105[) by cy e , then taking the sum over j = 
1, . . . , n and integrating the equality over (0, t), we have 

2||H^llL- (fl) +4jf i ^VxH*.«||} | , (fl) d. 

-'iff ((H fc ' e ) T Vu' £ ' e H fe ' £ + iv(|H fc ' £ | 2 )-u fe ' £ )da;)d S 
Jo JQ 2 

= 2||r»Ho||}, (114) 
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Summing up (|113[) - (lll4p . we obtain 

^\\b im + \\y^^\\b (a) +m^ 2 



L («) 1 v ' M L (Q) 11 "L {a) 

t 



4 / ((S^,D(u^)+H|VxH^||^ (fi) jd s 



'0 



<2||^r™uo||^. r +2||H ||^ (115) 



(n) 1 "«L (n) 
Using Korn's inequality to the fourth term we get 



< C*. (116) 

Multiplying the jth equation in (|104[) by &' e , then taking the sum over j = 
1, ...,n, using (|112l) with z — ^ „ I and integrating the equality over (0,i), 



we have 



< cwvfiiss \\b m + C £ ||S*- • D(u*.')|||, (fl) * 

||V x H^Uj^dt < C(n). (117) 

For fixed fc € N we can multiply the jth equation in (|103[) by , the jth 

equation in (|104p by and the jth equation in (|105l) by . Then we 

obtain 





da* 


r T db ke r T dc k ' e 




Jo 


dt 


- ?d t, I | ^ l-dt, I | ^ |^<C(n). 


(118) 


By miop-mnii, 


we 


can have the following convergence result as e — 


-> 






e — > p k * -weakly in L°°([0, T] x ft), 


(119) 


a fc ' e 




weakly in W 1,2 (0,T) and strongly in C([0, T]), 


(120) 


6 fe,e _, 


6 fe 


weakly in W 1,2 (0, T) and strongly in C([0, T]), 


(121) 




c k 


weakly in W 1,2 (0,T) and strongly in C([0,T]). 


(122) 



It follows from ([TT9"|) and (fT2"2"j) that 

u k,e > u k strongly in L 2r(Q 5 j,, W h2r^^ (^3) 
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jjfce > jjfc strong l y in L 2r(Q ; rp. W l>( fi )), ( 12 4) 

By (fTT9|) . (fl23)) - ([T24t . we can take the limit in the weak formulation of (fT02|) 
and get the transport equation 

p\ + div(p fc u fc ) = 0, on [0, T] x Q. 

Then applying Diperna-Lions theory (see [TU1I5U] ) of the renormalized solutions 
to the transport equation, we conclude that 

p k,e — y p k strong m L 2 ([0,T];L 2 (fi)) a.e. in (0, T) x Q. 

Thus by the above convergence results, we can take the limit e — > and 
obtain the solution (p k , u fc , 9 k , H fc ) solving the following system 

p\ + div(p fc u fc ) = 0, on [0,T] x Q, (125) 
(/^u fe ,^) + (p k [n k ]u k ,^) + (Sfc.Dtyj)) = ((V x H fe ) x H fe ,^), (126) 

(p k j t Q(6 k ),Lo 3 ) + {p k Q{9 k )xi k ,u. J ) + (q( Pk ,6 k , W fe ),Vw,)) 

= (|V x H fe | 2 ,^) + (S fe ,D(u fc H),(127) 

(—B. k ,vjj) + ^(curlH fe ,curlro i ) - (u fe x H fe ,curlro i ) = 0, (128) 

with 

p k -n = on [0,T], p„ < p k < p* , p k (0,.)= Po in (0,T) x Q, 

u fc (o,-) = r n u , 8 k (o,-) = r k e 0l H fe (o,-) = r fe H , 

where 

S k :=S(p k ,e k maxl T> k ) with e k max :=m*K{0\Q}. 



5.2 Limit k — > oo 



In this subsection, we establish some uniform estimates with respect to k. We 
follow the method in [7UT7]. Using the similar procedures of (fTTo) , ([TTT]) . (fTTg)) 
and (|125p , we have 

P*<p fc </A J WptW^^dtKC Vg 6 (1,00), (129) 



sup {||u fc ||^ 
te[o,T] ^ 



H 



VxH fc ||^ (f2) )rfi<C,(130) 
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Using Lions-Diperna theory of renormalized solutions, by (|129p . we have 

p k — >p strongly in C([0, T];L P ((2)) a.e. in (0, T) x Q. (133) 
By l(T5tj)l and 1(153)1 . we get 

u fc — >u strongly in L 2r (0, T; W*' 2r (J?)), 

and 

H fc ^H strongly in L 2 (0, T; W* ,2 (/2)). 
It follows from (fT^gj) - (fT^Tj) that 

sup {||(p^)(t)|| 2 }+ / |V0 fc | 2 «^ 
te[o,T] ^ 7q 



/ (6> fe ) a |Ve» fe | 2 da;ds < C(n), (134) 



0» 2 for 9 < 6*. 



+ 

/((0,T)xfi)/Q 

where Q = {(t,a;) G (0,T) x Q;0 k (t,x) < 6»„}. 
Define 

- / r for 6 > 6* 
1 ' ■ 1 0? /or < 0, 

and 

Using (|134[) . we have 

SU P H^L 2 ^ + / T |V^)| 2 ^ < C(n). (135) 
Using the similar procedure in [7] (also see Q7]), by (|135[) , we get 

^ sup 11^)11^ + V^n)* ^ C ( n )' ( 13fi ) 

/ \\ c F\\T' m rn\dt < C(n) with to = 2 /or a < and to = ^ + — /or a > 0, 
J ^ 3a + 5 
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\\(n(6 k )) 2 Kfc|lL 9 ((o,T)xf2) - G ( n ) with q = oo for a <0 

and q = 2(3a + 10) for a > 0, (137) 
3a 

/ \V9 k \ q dxdt < C(n) with q = 5(/ * + ^ for a < and q = 2 for a > 0. 
Jo a + 5 

Then it follows from the above estimates, (|127[) and the continuity of the 
projection r that 

\\*Wh< l0 ,TW* m * c V) f» ^t^- (138) 

By (|131D - (|138I) . we can establish the following convergence results as k — > oo 
9 k ^9 weakly in L«(0,T;W 1,9 '(«)), (139) 

p k 6 k ->> P 9 * -weakly in {z e L°°(0, T; L 2 (fi)), z t € l/(0,T; W _1 '''(fl))}. 
Consequently, using Aubin-Lions Lemma 1, we have 

p fe & — -> p9 strongly in C(0,T; (W 1,9 (/2))*). 

Then 

(pM> (W i.« (n)) . = I (p 9 ' 9 ^ 

J 

which together with (|139l) shows that 

9 k — >9 strongly in L 2 ((0,T) x Q) a.e. in (0, T) x Q, (140) 
This combining with ffl33|) . (fT36)) -([T37 | implies that 

K(9 k )^K(9) weakly in L 2 (0, T; W 1:2 (J?)), 

(R(e k ))-?K k — > (k(0))-3k strongly in L q " (0, T; L 9 * (.!?)) Vg*< ? , 
and hence 

q fc q := ft(p, 6»)V6" weafcZy in L*(0, T; W 1,(5 (J?)), 

where k := K(p,9 max ) with 6» max := {0, 0*}. 

Therefore, the above convergence results allow us to take the limit in (|125jl - 
(H"2gf and to obtain 



/ (p t , z) - (pu, Vz)dt = 0, p*<p<p* 
Jo 
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for any z € L 9 (0, T; W 1,S (J?)) for any q e [l,oo). 

((pu) t ,^-)-(pu®u,V^) + (S(p,0,D(u)),D(^)) 

+ / H T V^H + iv(|H| 2 ) • <^-cte = 0, 

/ (((pQ(0)) t ,h) - (pQ(6)u,Vh)-(d(p,9,V6),Vh))dt 
Jo 

= [ xH\ 2 ,h) + (s( P ,e,D(u)),n(u)h))dt, 

Jo 

for all h e L°°(0,T; W 1,? ) with g sufficiently large, 

T 

((H t , 6) + i/(curlH, curl6) + (u x H, curlfc)) dt = 0, 

for aU6eL 2 (0,T;W 1,2 (fl)). 
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